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^ ' Abstract 

o 

Q^ , We calculate the coefficients of the operator product expansion (OPE), in Polyakov's 

approach for Burgers turbulence. We show that the OPE has to be generelized and it 
Oh! 

(— I I is shown that the extra term gives us the instanton solution (shock solution) of Burgers 

•rH . equation. We consider the effect of the new-term in the OPE, on the right and left- 

X 

' tail of probability distribution function (PDF). It is shown that the left-tail of PDF, 

where is dominated by the well-separated shocks behaves as W{u) ~ n^'^/^. Finally we 
calculate the assymptotic behaviour of the N-point generating function of the velocity 
field, using the new OPE. 
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1- Introduction 

A theoretical understanding of turbulence has eluded physicists for a long time. A sta- 
tistical theory of turbulence has been put forward by Kolmogorov [1], and further developed 
by others [2-4]. The approach is to model turbulence using stochastic partial differential 
equations. In this direction, Polyakov [5] has recently offered a field theoretic method to 
derive the probability distribution or density of states in (l-l-l)-dimensions in the problem 
of randomly driven Burgers equation. The importance of the Burgers equation is that, it is 
the simplest equation that resembles the analytic structure of the Navier-Stokes equation, 
at least formally, within the scope of applicability of the Kolmogorov's arguments [6] . 

Polyakov formulates a new method to analyse the inertial range correlation functions 
based on two important ingredients in field theory and statistical physics, the operator 
product expansion (OPE) and anomalies. He argues that in the limit of high Reynolds 
number because of existence of singularities at the coinciding point, dissipation remains 
finite and all sublleading contributions vanish in the inertial range. Using the OPE one can 
find the leading singularities and show that this approach is self-consistent. 

Using the OPE Polyakov reduces the problem of computation of the veloctiy correlation 
functions to the solution of a certain partial differential equation. In this direction the N- 
point generating functions have been found in [7] . The two-point functions for different type 
of correlation for noise has been found in [8] and the perturbative calculation in the presence 
of pressure has been done in [9] . 

Here using the results of [7], we calculate the OPE coefficients proposed by Polyakov, 
and show that we have to modify the OPE. We write the generelized OPE and find the first 
correction to the PDF which was initially found by Polyakov. It is shown that the extra term 
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in the OPE gives us the instanton, which has found by Gurarie, Migdal [10] and Bouchaud, 
Mezard [11]. 

We find the right and left-tail of the PDF and show that the right-tail remaines unchanged 
and the left-tail , where is dominated by the shocks behaves as W{u) ~ which is in 

agreement with recent numerical experiments [12]. Finally we find the assyptotic behaviour 
of the N-point generating function with the new-term (instanton-term) in the OPE. 
2- The Instanton Solution of the Burgers Equation 

We consider the Burgers equation in one-dimensin: 

Ut + UU^ = uu^^ + f{x,t) (1) 

where u and u are velocity field and viscosity, respectively. The f{x, t) is a Gaussian random 
force with the following correlation: 

< fix, t)f{x\ t) >= k{x - x')5{t - t) (2) 

To investigate the statistical description of eq.(l) following Polyakov [5], we consider the 
following generating functional 

N 

^iv(Ai, A2, . . . Xn, xi,... xn) =< exp(^ Xju{xj, t)) > (3) 

Noting that the random force f{x,t) has a Gaussian distribution, satisfies a closed 
differential equation provided that the viscosity 1/ tends to zero: 

d 1 dZ 

+ I] ( y-Q^) = J2Hxi- Xj)XiXjZN + Dn (4) 

where Djv is : 

-Dat = i/^Aj < u" {xj,t)expY^Xku{xk,t) > (5) 



The first term in the r.h.s. of eq.(4) can be derived by applying the following relatin and 
then using the Gaussian character of the noise term in time and also by giving address to 
causality. 



However to remain in the inertial range we must keep v infinitesimal but non-zero. 
Polyakov argues that the anomaly mechanism implies that infinitesimal viscosity produces 
a finite effect. To compute this effect, Polyakov makes a conjecture for existence of an oper- 
ator product expansion or the fusion rules. The fusion rule is the statement concerning the 
behaviour of correlation functions, when some subset of points are put close together. 

Let us use the following notation; 



where ex^{xi) = e'^i"^^!). 

The Polyakov's conjecture is that in this case the OPE has the following form, 

d 

eM{x + y/2)ex^{x-y/2) = A(Ai, A2, y)eAi+A2(2;) + -B(Ai, A2, y)^eAi+A2 + o(y^) (8) 

This implies that Zj^ fuses into functions ^jv-i as we fuse a couple of points together. 
This conjecture allows us to evaluate the following anomaly operator (i.e. the Dj^-ierva. in 





Z(Ai, A2, . . . . . .xjv) =< eAi(xi) . ..ex^{xN) > 



(7) 



eq.(4)). 




(9) 



which can be written as: 




(10) 
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As discussed in [5] the possible Galilean invariant expression is: 

ax(x) = a{X)ex{x) + P{X)^ex{x) (11) 

Therefore in steady state the master equation takes the following form, 

d d ~ 

m = m + \ (12) 

Let us consider following correlation for / in /c-space as follows [13]: 

< f{K t)f{k',t) >= ^ko5{e - ^)5{k + h!)5{t - 1) (13) 

therefore we obtain: In the inertial range where Xi — Xj « L, we find: 

k{x, - X,) = ko{l - ^^^^^) (14) 

Polyakov has found the following explicit form of Z2 for k{xi — Xj) given by eq.(14): 

Z2{liy) = el^>^y^"' (15) 

and the following expression for density of states as the Laplace transform of ; 

W{u,y)= ^e-^^Z^ii^y) (16) 

where = 2(Ai — A2) and y — x\ — X2- Now one can write the right and left-tail of W{u^ y) 
as following: 



e ^^y' if {u/y) — > +00 

W{u,y) = { (17) 
^3/2^-5/2 ^ ^9/2^-11/2 if ^ _^ 



In ref.[7] we have found the exact N-point generating function (i.e. the Z^) as follows: 

Zn = (A1A2 • • • Ajv)''^(At2At3 • • • /x7v)"^^e'/=^('^^^^+^3f3+-+'^~^^)'^' (18) 
where — ^2N^ ^» given by: 



X1+X2 + X3 + ...Xn 

yi = 



Z/2 = Xi 

ys = X2- 



N 

X2+X3 + ...Xn 

A^- 1 

X3 + X4 + . . .Xn 



N-2 

and un — xn-i — xn (19) 



and 



Ai + A2 + • • • + Ajv 
= 

A^-1.^ A2 + A3 + ... + A;vi 
= N^l 1 

N-2 A3 + A4 + ... + AJV1 

and /In — 2(Ajv-i — Xn) (20) 



It follows that the A^-point correlation function of v is: 

N 



GW(xi, • • • , x^) - Zim.^oA-^ E (21) 

k=0 



where a^^^ are some constants. 

Now we try to find the OPE coefficients in eq.(8). To do this according [7] we use Z^ and 
tend X2 close to X3. Z2 is given by eqs.(15) with Ai + A2 = 0, and Z3 has following form: 

Zs = (AiA2A3)^/6[3(Ai-^^^)(A2-A3)]-^/^ 

^{2/3[3/2(Ai - - f2+f3 )+2(A2 -\3)(x2 -xa)]^/^} 



(22) 
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where Ai + A2 + A3 = 0. In eq.(22) we take — X2 — 2e and it is easy to show that in 
the hmit e — > we find: 



Z3 = (AiA2A3)5/6(9/2(A^ - A2))-5/4e{2/3[2(Ai-{A.+A3)){.i-x.)]3/2} 

{ ! + (-! + 7/64 + • • •)7/8(A2 + A3)(-2)[-4(A2 + \^)Y'\x, - x^f/^ 

+ (1 - 7/32 + •••)^^^^(-2)[-4(A2 + A3r^(x,-X2)^/^6 + 0(6^)} (23) 

Comparing with eq.(8) and using eq.(15) we find that: 

A{X2, A3, e) = (-A2A3(A2 + A3))^/^(^(A^ - A^))"! (24) 

7A — 9'^A 

S(A2, A3, 6) = (1 - 7/32 + . . •)^^^^^^^(A2A3, 6)6 (25) 

Indeed we determine the A(A2, A3, e) and B{X2, A3, 6) by means of the first and third terms 
in eq.(23). However the second term in eq.(23) (i.e. which is porpotional to {xi — 0:2)^^^), 
can be treated as ^-^2 and this will change the OPE (i.e. eq.(8)) to the following form: 



d 

ex^ix + e/2)exs{x - e/2) = ^(A2, A3, 6)eA2+A3(2;) + 5(A2, A3, 6)— eA^+Ag 

+ C(A2,A3,6)^eA,+A3+0(6^) (26) 

where using eq.(54), C(A2, A3, e) can be written in terms of ^(A2, A3, e) as follows: 

C(A2, A3, 6) = (-1 + 7/64 + . . .)I(A2 + A3)A(A2, A3, e) (27) 

Using eq.(26) one can show that the modified OPE leads to following equations for ax{x) 
and Z2 (i.e. eqs.(ll) and (15)): 

ax{x) = a(A)eA(x) + ~P{X)—ex{x) + 7(A)— eA(x-) + • • • (28) 



{d^ - —)dyZ2 + cnd^Z2 - plVZ2 = (29) 
A* 

where we have used /9(A) = ^ , a — according [5]. We have used the scahng arguments 
and show that 7(A) has the following form: 

7(A) = cA (30) 

Therefore we find following assymptotic behavior for Z2: 

Z2{iiy) = ei(^^)'^'-t('^^) (31) 

First of all it appears that in the hmit c — > 0, one finds the Polyakov's result for Z2. 
Also by using of eq.(28), one can show that the 1-point function of the velocity behaves as 
< v{y) >~ y, which is related to the well-known shock wave structures [6], and is exactly 
the instanton solution of the noisy Burgers equation, found by Gurarie and Migdal [10], 
and later in a different way by Bouchaud and Mezard [11]. The interesting point is that 
Polyakov's approach generates in self-consistent way the instanton solution. 

Let us write the right and the left-tail of the PDF. It is easy to show that the W{u,y) 
behaves as follows: 



W{u,y) = < 



e 3S"^2^ if (u/y) +00 

(32) 

^3/2^-5/2 ^ cy^/2^-y^ + c2y7/2^,-9/2 ^ ^9/2^-11/2 ^ _^ 



It is evident that the right-tail remanies somewhat unchanged and the left-tail of the 
PDF, where is dominated by the shocks behaves as W{u) ~ which is in agreement 

with the observations [12]. 
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Now let us find the assymptotic behaviour of the N-point generating functions. For 
investigating the assymptotic behaviour we use the following important property that: 

^l.Xid^, = (33) 

where /ii — 0, which in turn gives the following assymptotic behaviour of the N-point 
generating functions in the limit of E^j^/^jt/j =^ oo: 



Now it is easy to show that for un- forced Burgers equation {ko — 0), we find, 

< w«(r) >~ r« (35) 
which is again cosistent with other approaches [14]. 

Acknowledgements: I would hke to thank A. Rastegar, S. Rouhani and J. Davoudi for 
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